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Combination of the Liouville equation with the g-averaged en- 
ergy U q — {H) q leads to a microscopic framework for nonex- 
tensive ^-thermodynamics. The resulting von Neumann equa- 
tion is nonlinear: ip = [H,p q ]. In spite of its nonlinearity 
the dynamics is consistent with linear quantum mechanics of 
pure states. The free energy F q — U q — TS q is a stability 
function for the dynamics. This implies that g-equilibrium 
states are dynamically stable. The (microscopic) evolution 
of p is reversible for any q, but for <; / 1 the corresponding 
macroscopic dynamics is irreversible. 

Standard thermodynamics is based on the Gibbs- 
Shannon-von Neumann entropy Si = -fc^Tr (plnp) and 
the internal energy U\ = Tr pH. The equilibrium den- 
sity matrix po minimizes the free energy F\ = U\ — TS\ . 
A microscopic foundation for thermodynamics is based 
on the von Neumann equation (vNE) 



ip = [H,p\. 



(1) 



Equilibrium states are stable fixpoints of this dynamics. 

It is perhaps not so widely known that the vNE can be 
regarded as a classical Hamiltonian system with Hamil- 
tonian function U\ . In this context F\ is a stability func- 
tion for the underlying Hamiltonian Lie-Poisson dynam- 
ics. We shall elaborate on these points later but first we 
want to pose the following problem. It is known that 
there exist physical systems that are naturally described 
by a nonextensive thermodynamics |lj . The structure of 
this theory is analogous to the ordinary one, with one ex- 
ception: Instead of U\, Si and F\ one takes their q 1 
generalizations U q , S q and F q (see below). The ques- 
tion is whether there exists an underlying dynamics for 
nonextensive thermodynamics. Is it given by the stan- 
dard vNE (as worked out in [§)) or, maybe the dynamics 
should also be g-modified? 

The analysis given below is based on one single Ansatz: 
the g-averaged energy U q is also the Hamiltonian func- 
tion of the quantum system. As we shall see, this im- 
plies that F q is again a stability function for the micro- 
scopic dynamics. In particular, the equilibrium states of 
^-thermodynamics are dynamically stable. 

The vNE is an immediate consequence of the 
Schrodinger equation (SE) 



(2) 



if p represents a pure state, i.e. p = g = \vj) (we use g 
to denote pure state density matrices or their reductions 



to subsystems). However, in real experimental situations 
one does not deal with pure states since (a) there exists 
a classical lack of knowledge about quantum sources and 
(b) entangled states lead to non-pure density matrices 
when reduced to a subsystem. Indeed, if the classical 
state of the device that prepares the quantum ensem- 
ble is not exactly known, and the experiment is repeated 
several times, one has to apply a purely classical averag- 
ing over the classical configurations of the source. This 
results in a non-pure density matrix 



(3) 



p= / dgw{g)g 



where one integrates over quantum-mechanical pure 
states whose distribution is given by w(g). p can be used 
to calculate experimental averages since the latter are 
linear in g: 



Tr pA ~ J dg w(g) Tr gA. 



(4) 



The distinction between the state g and the distribu- 
tion w(g) (function defined on the space of states) leads 
ultimately to a sharp distinction between the Liouville 
equation (LE) and the vNE. In order to explain this main 
point let us first note that the vNE ([!]) can be written as 
a classical Hamiltonian system || [?) : 



ipa = {p a , (H)x} 



(5) 



where p a = are the components of p taken in 

some basis, (H)i = Tr pH is the average energy, and the 
bracket 



{ A,B} = Pa n\ c ^^-. 

opb opc 



(G) 



is a Lie-Poisson bracket on the manifold of states ||. 
Here il a bc are structure constants of a Lie algebra 
[gl(n, C) for finite dimensional Hilbcrt spaces] and the 
summation convention means summation or integration 
with respect to an appropriate measure PJlCfl . Using the 
standard argument one derives the LE 



iw = {w, 



(7) 



It is obvious that the linearity in w of the LE (J7|) is com- 
pletely unrelated to the linearity of the vNE, and follows 
directly from the fact that the dynamics (||) is Hamilto- 
nian. Also the physical meaning of the linearity of (0) is 
clear: It reflects the linearity of averaging and the fact 
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that an experimentalist can control the form of w by im- 
proving the measurement device. 

The second class of density matrices that are non-pure, 
in the sense that g 2 ^ g, occurs as an entirely quantum 
phenomenon and is a result of entanglement between cor- 
related quantum systems. Such states are fundamentally 
and irreducibly mixed. They can be written in different 
ways as convex combinations of pure states and all such 
combinations have to be regarded as physically equiva- 
lent HUH]. 

On the other hand, the decomposition (^) is uniquely 
determined by the experimental setup. The difference in 
physical status of the "pure-state" decompositions of p 
and g implies that there exists a physical difference be- 
tween the linearity of the LE (0) and this of the vNE (0) . 
The linearity of the latter is a postulate that is indepen- 
dent of both the linearity of (Q) and this of the pure state 
SE (||) . Linear SE is compatible with any equation of the 
form 



ig=[H,f(g)} 



(8) 



provided f(g) = g for g 2 = g, which holds for all func- 
tions satisfying /(0) = and /(l) = 1. The choice of 
f(x) — x is convenient but does not seem to be dictated 
by any fundamental principle. On the contrary, we will 
argue that other choices of / may be physically relevant 
and, in particular, we will show that there exists a link 
of fix) — x q with the nonextensive g-statistics intro- 
duced by Tsallis Q [note that q > is needed to ensure 
/(0) = 0]. 

We shall first show that the modified dynamics given 
by Eq. (|J) has the same Lie-Poisson structure as ([!]). 
Assume / has a Taylor expansion f(x) = Yl'kLi fkx k 
with radius of convergence at least 1. Consider the Lie- 
Poisson dynamics 

ig a = {ga,(H) f } (9) 

with the 1-homogeneous Hamiltonian function 



(H) f = Tr (Trp)/ 



Q 



Tig 



H 



(10) 



Variation of (H) f with respect to g gives the effective 
Hamiltonian 



H(g) = 6 Ml = ^MTrg)^ £ g^^Hg- 



fe-i 



5q 



k=l 



m— 



where /' = df/dx. (^) when written in an operator no- 
tation is 



ig=[H(g),g] = (Trg)[H,f(g/Trg)]. 



(11) 



Tr g is an integral of motion so we can consider solutions 
normalized by Tr g = 1 which shows that M) is indeed 



a particular case of ([)]) . 1-homogeneity of (H) j implies 
that (H) f = Tr gH(g). Taking f(x) = x q , a normalized 
g, and denoting the corresponding (H) / by (H) q , we can 
see that (H) q = Tr gH{g) = Tr g q H i.e. the average ef- 
fective energy equals the g-average of H, an internal en- 
ergy typical of Tsallis generalized thermodynamics. Such 
averages were shown to be naturally linked to nonexten- 
sive g-entropies 

ff-H. Obviously, any theory that for 
some reasons deals with q-averages involves some degree 
of nonextensivity independently of what kind of entropy 
is used. Let us note that in the above case the vNE is 
nonlinear, 



ig=[H(g),g}= [H,g q ]. 



(12) 



Nonlinear vNE has been studied recently in |20| , and is 
used regularly in statistical physics. The nonlinearity is 
usually due to friction forces and should be compensated 
by adding a noise term to the vNE in order to keep the 
average energy constant. However, for a nonlinearity of 
the form (^) the energy is a conserved quantity. Hence 
there is no need for the balancing noise term. vNE's 
of the form @ were independently found in the con- 
text of a Lie-Nambu dynamics and studied in |9|j2l||. It 
was shown, in particular, that their Hermitean Hilbert- 
Schmidt solutions possess time-independent spectra, an 
important fact that allows to treat the solutions g(t) 
as density matrices. This implies also that for any t\ 
and t2 there exists a unitary transformation satisfying 
e(*a) = U{t 2 M)Q{ti)U{t 2 M)~ l - Although U(t 2 ,h) ? 
U (1,2 — ti, 0) [equality would imply linearity of evolution] 
the local generator of U(t + e, t), for e — > 0, exists and is 
our effective time-dependent Hamiltonian H[g(t)). De- 
note C n = Tr(g n ), n e N. One finds {C n ,F} = 
for any F, which shows that C n are Casimir invariants 
for the dynamics. The set of invariants contains also all 
n-averages of H since {(H) n , (H) m } = for any nat- 
ural n and m. The dynamics we consider is therefore 
so regular and so close to the linear one that one may 
wonder whether such equations do possess nontrivial so- 
lutions. Fortunately, the answer is positive. An analyt- 
ical Darboux-type technique of solving (12) for q = 2 
has been recently developed Q and various explicit so- 
lutions were found. In a classical context the q = 2 case 
was discussed in great detail in relation to Lie-algebraic 
generalizations of classical Euler equations (cf. || and 
references therein, in particular [|23|,^4|). It is also well 
known that similar q — 2 Lie-Poisson equations describe 
plasma dynamics (in the context of a generalized statis- 
tics a paper of particular relevance is P5| ) . 

Since the description we propose is meant to provide 
a fundamental quantum background for a generalized 
statistics it must be also capable of dealing with collec- 
tions of nonextensive systems. This means we have to 
provide a recipe for extending the von Neumann dynam- 
ics from subsystems to composite systems. The exten- 
sion should be self-consistent in the sense that a dynam- 
ics of a subsystem should be independent of whether the 
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system is considered alone or as part of a collection of 
many noninteracting systems. This is achieved by taking 
the 2-system Hamiltonian function {Hi) 1 + {Hn) q ^. The 
(^-averages occuring in the Hamiltonian function of the 
composite system are 

{HiY qi = Tn({QirHi), 
{Hii)ll = Tr ii{{qii)^Hu) 

and the density matrices gi and qh are the reduced den- 
sity matrices of the respective subsystems. The two- 
system equation one obtains is 

iQi+n = [Hi{gi) ® 1// + 1/ ® Hii(qh), qi+ii], (13) 

where Hi{gi), Hn{gn) are the effective Hamiltonians 
of the subsystems. The choice of this particular type 
of extension follows from the general results proved for 
a Lie-Poisson nonlinear quantum mechanics of density 
matrices in ||| 10 2fJ. The expressions C n {gi + n) = 
Tr i + ii{gi + u) n are time-independent (as Casimir invari- 
ants) for any natural n. Therefore, if gi+n{i) is a Her- 
mitean Hilbert-Schmidt solution of (|13| ) then its eigen- 
values are time independent on the basis of the standard 
argument |pl| , and therefore gi+n {t) is a density matrix 
if it is one at t = 0. On the other hand, taking partial 
traces of (|l3| ) one verifies that 

igi = [Hi, {gi)^}, ig H = [H Ul {gn) q % 

as required by the self-consistency of the extension. All 
these results have an immediate extension to more gen- 
eral Hamiltonian functions {H)f. 

Having established all these general results we are now 
in position to discuss in more detail the links to the gen- 
eralized statistics proposed by Tsallis 0. It is based on 
the internal energy U q = Tr g q H and the corresponding 
entropy 



S q {g) 



Tr g — Tr {g q 



(14) 



U q is naturally associated with S q since then various re- 
lations typical of q = 1 thermodynamics turn out to be 
q- independent. However, standard thermodynamics is 
static and the relations between U q and S q are evalu- 
ated in thermal equilibrium. From the dynamical point 
of view an equilibrium state go is a fixed point of the dy- 
namics, i.e. [H, go] — 0. There exists an infinite number 
of such states but not all of them have to be stable if a 
nonlinear Lie-Poisson dynamics is involved. The stabil- 
ity tests that are typically used in such a situation (say, 
in hydrodynamics and plasma physics) are the energy- 
momentum, energy- Casimir [ p7[ or energy-invariants pq ] 
methods. In the energy- Casimir method (used when one 
knows the Casimirs but does not control the symmetries) 
one looks for minima or maxima of the "stability func- 
tion" 



where h{g) is a Hamiltonian function of the Hamiltonian 
dynamical system and $ is a function of the Casimir 
invariants Ck typical of this system. The latter function 
is determined by the requirement that X{g) has a strict 
minimum or maximum at go, in particular 



SX , s 
^(,o)=0. 



(16) 



In our case h{g) is Uf = {H) / (or U q if we restrict the 
analysis to f{x) = x q ), and the Casimirs are all functions 
C{g) that can be written as a trace of a convergent power 
series i.e. 



C{g) = Tr 



c fe C fe (0) -^(d, C 2) ...). 



It is clear that the stability function X for the energy- 
Casimir method is nothing else but the free energy F 
corresponding to a generalized entropy S = — $/T with T 
the temperature. In this way the thermodynamic relation 
F = Uf — TS is recovered. The equilibrium state go is an 
extremum of F. If this is a strict minimum (or maximum) 
then the orbits of density matrices in a neighborhood of 
go are dynamically stable. Thermodynamic stability of 
the Tsallis thermodynamics has been raised in (29) and 
settled in poj . Therefore equilibrium states extremizing 
F will generically be dynamically stable fixed points of 
the nonlinear vNE. 

Once dynamic stability of go is established it becomes 
meaningful to study linear response theory pl| . This has 
been done in the context of non-extensive statistics by 
Rajagopal 17 3^]. However, the theory has to be modi- 
fied because of the nonlinearity of the vNE. A discussion 
of these modifications is out of the scope of the present 
paper and will be presented elsewhere. 

Let us illustrate the above results with the simple ex- 
ample of a single spin in an external field. The Hamil- 
tonian is given by H = —j-icr z (assume fi > 0, the 
a a , a = x,y,z are the Pauli matrices). A general Her- 
mitean 2x2 matrix with eigenvalues Ai, A2 is 

q = -(Ai + A 2 )l + \{X X - A 2 ) cos(</>)er z 



(Ai - A 2 ) sin(0) {cos(4>)a x + sin{ip)a y ) . (17) 



If sin(0) = or Ai = A 2 then g is time- invariant. In the 
other case Eq. (0) implies = 0, 



-lu with lo — 2/i 



/(Ai)-/(A 2 ) 
Ai — A-2 



(18) 



X(Q) = h(Q) + $(d,C 2 ,...) 



(15) 



and the Larmor precession frequency to depends in gen- 
eral on the eigenvalues of g. For this example this is the 
only effect induced by the nonlinearity of the vNE. Note 
that lo = 2/i is still valid for pure states — the underlying 
quantum mechanics remains unchanged. 

One verifies that {H) q = -fi (A« - (1 - A) 9 ) and 
H{go) — — + a multiple of 1 with 7 = fJ.q[\ q + (1 — 
A) 9 l/2. Now assume f{x) — x q . The internal energy is 
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U q = Tr g q H = —fx cos(^) (A? 



A!) 



The g-entropy equals 5 9 = fcs [l — A^ — A|] /(?— 1) In both 
formulas one should take Ai = A > and A2 = 1 — A > 0. 
Of physical interest are minima of F for which A > 1 — A 
and U q < 0. Hence one can take cos(</>) = 1. From 
dF/dX — 0, with j3 = l/fc^T and assuming < \q — 
l\/3fi < 1 one finds that the eigenvalue A of g is the 
solution of 



A 



1 - A 
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The value of d 2 F/d\ 2 at equilibrium is strictly positive. 
In fact, F has an absolute minimum at g = go. This 
thermodynamic stability implies dynamic stability of go 
as a fixed point of the nonlinear vNE. One concludes that 
Tsallis thermodynamics is useful to analyse the dynamic 
stability of fixed points of the nonlinear vNE. 

A special feature of the nonlinear vNE is that classi- 
cal mixtures of initial conditions evolve irreversibly with 
time (this property has been discussed extensively in 
@). Take e.g. 



1 />7T /■2-7T 

dA / sin(0) dcj) / dip iu(A, <j), i)/)q{X, <j>, ip) 
Jo Jo 



with w(X, 4>, = (1/8) sin(^/2) and with g(X,<j),ip) 
given by (|l7|). A short calculation then shows that 



Pit) 



24 



dA(2A - 1) 



x (^cos(uj(X)t)cT x + sin(w(A)i)er, y 



with w(A) given by (jl8|). Due to the dependence of uo 
on A a dephasing occurs and the classical density matrix 
converges to |l as t — * 00. The lack of knowledge about 
initial conditions leads to a true irreversible decay. 

An analysis of more complicated examples including 
linear response theory will be presented in a forthcoming 
paper. 
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